Quantum dimer model with Z 2 liquid ground-state: interpolation between cylinder 
and disk topologies and toy model for a topological quantum-bit 
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We consider a quantum dimer model (QDM) on the kagome lattice which was introduced recently 
[Phys. Rev. Lett. 89, 137202 (2002)]. It realizes a 7,2 liquid phase and its spectrum was obtained 
exactly. It displays a topological degeneracy when the lattice has a non-trivial geometry (cylinder, 
torus, etc). We discuss and solve two extensions of the model where perturbations along lines 
are introduced: first the introduction of a potential energy term repelling (or attracting) the dimers 
along a line is added, second a perturbation allowing to create, move or destroy monomers. For each 
of these perturbations we show that there exists a critical value above which, in the thermodynamic 
limit, the degeneracy of the ground-state is lifted from 2 (on a cylinder) to 1. In both cases the 
exact value of the gap between the first two levels is obtained by a mapping to an Ising chain in 
transverse field. This model provides an example of solvable Hamiltonian for a topological quantum 
bit where the two perturbations act as a diagonal and a transverse operator in the two-dimensional 
subspace. We discuss how crossing the transitions may be used in the manipulation of the quantum 
bit to optimize simultaneously the frequency of operation and the losses due to decoherence. 



(N 
> 

m 

OS 

SO 
O 

o 



■ 

a : 
o . 
o ■ 



X 



I. INTRODUCTION 

Quantum dimer models (QDM)§j§ provide simple 
examples!™ of microscopic Hamiltonians with short- 
ranged resonating valence-bond ground-states (or dimer 
liquid) with gapped excitations and no broken symmetry 
(Z2 liquids). It has been known for a long time. that such 
liquids are characterized by topological orderH although 
the system breaks no symmetry and has no local order 
parameter, the ground-state acquires a degeneracy (in 
the thermodynamic limit) which depends on the genus 
of the surface where the model is defined (disk, cylinder, 
torus, etc.). Remarkably, such a topological degeneracy 
is insensitive to small Jojcal perturbations (such as weak 
disorder for instance) Blh3 On the other hand, it is clear 
that strong enough local perturbations should lift this 
degeneracy. Consider for instance a Z2 dimer liquid on 
a cylinder, with a two-fold degenerate ground-state. We 
turn on an external potential which penalizes (with an 
energy A > 0) any dimer sitting across a line extend- 
ing from one edge of the cylinder to the other. For very 
large A, one effectively "cuts" the cylinder down to a 
disk topology and one expects a single non-degenerate 
ground-state. It is therefore natural to expect a phase 
transition at some intermediate value A. 

We provide here a simple model where the spectrum, 
and the ground-state degeneracy in particular, can be ex- 
actly calculated as a function of A and the system size. 
This model generalizes a QDM on the kagome lattice 



(network of corner-sharing triangles with triangular and 
hexagonal plaquettes, see Fig. g)) which was introduced 
recently!] The full spectrum (and wave- functions) can be 
obtained in an elementary way and excitations consist 
of static and non-interacting Ising vorticesQ (visonsEJ). 
In this paper we show how the solution of the model ex- 
tends to a situation where an external potential is applied 
along a line of the system. The solution is obtained by 
noting that the bulk of the system decouples from the 
line and the line is exactly described by an Ising chain 
in transverse field (ICTF). As a result, we find a finite 
critical value A c of the perturbation below which the sys- 
tem behaves as a cylinder (the energy difference between 
the two quasi-degenerate ground-states is exponentially 
small in the system size). For A > A c the system behaves 
as a disk and the ground-state is separated from the first 
excited state by a finite gap 0(X — A c ). 

It has been argued that gapped systems with a topolog- 
ical degeneracy could provide physical ways to implement 
quantum-bits (qubits) which would be-pmtpcted from 
decoherence by their topological nature flDEJO Since no 
local measurement can distinguish the different ground- 
states if the system is infinitely large, any manipulation 
(unitary rotation) or measurement (projection) of the 
state of the qubit through local observables will have to 
rely on finite-size effects. We discuss this issue in Sec- 
tion [v| in the light of the present solvable model. The 
effective Hamiltonian acting on the two lowest levels is 
expressed in terms of two generators T x and T z of ro- 
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FIG. 1: A dimer covering on the kagome lattice (fat bonds). 
The corresponding representation with arrows living on the 
bonds of the hexagonal lattice (dashed lines) is displayed. 



tations of the qubit about two quantization axis. We fi- 
nally explain how one could take advantage of the phase 
transition at A = A c to perform unitary rotations. The 
problems of this approach, such as thermal excitations, 
will also be discussed briefly. 



II. ARROW REPRESENTATION, SOLVABLE 
QDM AND TOPOLOGICAL DEGENERACY 



• a h = l~L=i T i ■ Flips the 6 arrows i = 1 • • ■ 6-around 
the hexagon h (smallest closed loop on H ) .E3 From 
the definition of t z above, the a x operators satisfy 
(cr x ) 2 = 1 and commute with each other. 

• rf : Compares the arrow at site i with the the ar- 
row in some (arbitrary) reference covering (= +1 if 
they are the same, —1 otherwise). The hard-core 
constraint on dimers translates into Tjf rf rf = 1 for 
every triangle (012) of the kagome lattice. 

From now on and for most purposes one can forget the 
dimers themselves and. focus only on the bond degrees 
of freedom rf = ±l.lEl We note that, in principle, the 
rf = ±1 degrees of freedom could be physically realized 
with real spins living on a kagome lattice. A strong easy- 
axis anisotropy could then force them to point toward 
the center of one of the neighboring triangle. 



Bulk Hamiltonian 



The QDM introduced in Ref. | is 



wo = -E°s = -EII 7 '' 



hi 



(1) 



We consider a QDM defined on a kagome lattice with 
periodic boundary conditions in one direction (cylinder) 
but the arguments are easily generalized to other topolo- 
gies. 



All the eigenstates are easily obtained because the 
operators comrrutte with each other and have two eigen- 



values ai 



±iB 



A. Arrow representation 

Because it is the natural formalism to describe and 
solve the QDM discussed here, we begin by reminding the 
representation oLdimer coverings of the kagome lattice in 
terms of arrows.tSEj 

The sites of a kagome lattice K (noted i) can be iden- 
tified with bonds of the hexagonal lattice HS3 The tri- 
angles of K (noted t) are sites of H . As for hexagons of 
K (noted h), they also correspond to hexagons of H. 

From a (fully-packed) dimer covering of K we orientate 
the bonds of H (arrows) in the following way: Each bond 
of H is a site of K which has one dimer, the corresponding 
arrow points toward the interior of the triangle (of K) 
where the other end of the dimer is. This is illustrated 
in Fig. [l| As a consequence, the number of incoming 
arrow(s) is even (0 or 2) at each vertex of H . Inversely, 
any arrow configuration satisfying the parity constraint 
at each vertex defines a unique dimer covering. 

We can now define the operators t x , t z and a x act ing 
on the arrows. The notations are those of Refs. [| and 

• rf: flips the arrow at site i £ K. Any product 
T i T j ' ' ' aroun d a close loop (of H ) is a "physical" 
operator in the sense that it conserves all the con- 
straints. 



C. Topological sectors 

As usual for dimer models, the configurations are 
grouped in topological sectors (TS): two configurations 
are in the same TS if and only if they can he trans- 
formed into each other by a succession of loca£3 moves. 
As explained below, there are two TS when the system 
has the topology of a cylinder. 

First draw a cut A (crossing the bonds of the lattice) 
going from one edge of the cylinder to the other (Fig. ||). 
Let be the number of dimers crossing A. It has a 
simple expression in terms of the rf : 



"a = ±£(i- 



(2) 



where the sites i = • • • L are the centers of the bonds of 
H which are cut by A, as shown in Fig. ||. For simplic- 
ity, we assumed in Eq. ^ that no dimer crosses A in the 
reference covering. 

Any local dimer move conserves the parity of A^a and 
it is natural to define: 



T' t: 



n 



(-1) 



(3) 



All coverings with T x — 1 (resp. —1) define a TS, called 
S+ (resp. S-) and TIq can be diagonalized separately in 
each sector. 
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FIG. 2: Kagome lattice on a cylinder with a cut A going 
from one edge of the cylinder to the other. The dual cut A* 
passes through the centers of the triangles and winds around 
the cylinder. 




FIG. 3: Dimer covering of the kagome lattice in the vicin- 
ity of the cut A (dashed line). The arrows next to A (ap- 
pearing in Eq. ^) are shown. The signs in the hexagons 
Int.,/ii, • • • , /i_L,Ext. indicate the values of the correspond- 
ing pseudospins <r z (with the assumption that the reference 
configuration has no dimer crossing A). 



D. Topological degeneracy of Ho 

It is straightforward to check that Ho has the same 
energy in each sector. Let A* be a closed loop encircling 
the cylinder (Fig. ||) and define an operator 

T z =Y[t* (4) 

iSA* 

which flips all the corresponding arrows. T z commutes 
with all the a x operators and maps S+ onto 5_ : 

rpZrpX rpXrpZ /g\ 

This shows that if is an eigenstate of Ho, T z \ip) is 
an eigenstate of Ho with the same energy (but in the 
other TS). This demonstrates the two- fold (topological) 
degeneracy of the eigenstates of Ho- 



III. PERTURBATION LIFTING THE 
DEGENERACY BETWEEN THE T x = ±1 
SECTORS 

We introduce a potential energy term which couples to 
the dimers crossing A:EJ 

Hi = 2XN A (6) 

L 

= A5>"7f) (7) 

z=0 

As discussed in the introduction, we expect that a small 
A should not affect the two-fold degeneracy while A ^> 1 
should leave a single ground-state. In presence of H = 
Ho + Hi, T x = (— 1) N& is still a conserved quantity but 
T z does not commute with Hi and the two sectors are 
no longer degenerate when A ^ 0. Since for A ^ 1 the 
system minimizes Na, the ground-state of the T x = +1 
sector tends to a state with Na — and that of the 



T x = — 1 sector to a state with Na — 1. Instead of hav- 
ing a superposition of dimer configurations with different 
values of Na but a fixed parity (non-local observable), the 
large A limit corresponds to a well defined Na (sum of lo- 
cal operators). While T x is still a conserved quantity, we 
already see that the topological nature of the T x = +1 
and T x = — 1 sectors disappears when A is large. 

The perturbation Hi is identical to the one introduced 
by Ioffe et alQ in a triangular-lattice QDM in order to 
manipulate ("phase shifter") their qubit. However, in our 
case, the existence of an arrow representation makes it 
possible to calculate exactly the spectrum of H = Hq+Hi- 

A. Mapping to the ICTF 

The Hamiltonian H = Ho+Hi (Eqs. and 0) can be 
separated into a "bulk" and a "chain" part as below : 

H = Ha + Titmik (8) 

L L 

Ha = -E<+ A E(!- T n ( 9 ) 

i=l i=0 
Hbulk = - E a h' (10) 

It is important to emphasize that Ha and 7ibuik commute 
with each other and can therefore be treated separately. 
From now on we concentrate on Ha- 

a z pseudo-spins. — A o~\ operator can be introduced 
for each hexagon h in the following way. Due to the local 
constraint {t x t x t'£ = 1 on each triangle of K), the bonds 
of H where sf = —1 necessarily form non- intersecting 
closed loops.c3 We interpret these loops as domain walls 
for some Ising pseudo-spins er^ — ±1 which leave on each 
hexagon. To remove the two-fold ambiguity we assume 
that the exterior has a fixed spin a z Kt = 1. In turn, 
this defines a of nt associated to the "interior" (Fig. ||). 
It is easy to check that this Ising spin labels the TS of 
the configuration because of nt = JlieA* T f = T x ■ The 
other "bulk" pseudo-spins are those introduced by Zeng 
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FIG. 4: Cylinder (full lines) and loops (dashed lines) along 
which a dimer configuration c differs from the reference. The 
signs indicate the value of <r z in each domain. 



and Elserlll Eventually we mention that cr^ and afl an- 
ticommute (they commute if not on the same hexagon), 
as suggested by the Pauli matrix notation. This is eas- 
ily checked from the definition of cr^ in terms oX arrows. 
From the definition of ai we have the relation D'lj 



(see Ref. |l2] and Appendix [A]). On the other hand, in 
the ferromagnetic phase (A > 1), the Ising spins want to 
order and the boundary condition <7q ^ erf generates a 
finite energy cost AE ~ O(L ) (sec Eq. AT3). 

The result is thus that for A < 1 the ground-state 
is asymptotically two-fold degenerate in the thermody- 
namic limit. The gap between the two TS is AE ~ 
cxp(— L/£(A)) where £ is the correlation length of the 
ICTF. In the thermodynamic limit there is a finite criti- 
cal value A c = 1 above which the topological degeneracy 
is destroyed. Above A c the Ising pseudo-spins have a 
positive magnetization. Sincer-fr^ is an operator which 
creates an Ising vortex (visonla), it is natural to inter- 
pret (<r?) > as the existence of a condensate of those 
particles (along A). This condensation is responsible for 
changing the "effective" topology of the system from a 
cylinder into a disk. In this simple model what happens 
along the chain A is decoupled from the bulk of the sys- 
tem. The perturbation caused by the potential A does 
not extend into the bulk, which remains a liquid with 
non-interacting and static visons excitations. 



(11) 



where h and h' are the two hexagons touching i. For an 
arrow ig next to a boundary (say interior), this relation 
is modified to: 



These relations are used to get: 



Ti-A 



A £<< 

i=0 



(12) 



+ (L+1)A (13) 



where the two boundary spins are identified to cTq = erf nt 
and &h L+1 = o'oxt- 1^ appears that Ha is nothing but the 
Hamiltonian of an open ICTF with a magnetic exchange 
A and a transverse field equal to 1. This model can be 
solved by a standard Jordan- Wigner transformation and 
maps onto free fermions. In the thermodynamic limit it 
has a paramagnetic phase with (0%) = for |A| < A c = 1 
and an ordered phase with (a^) ^ for |A| > A c = 1. 

The boundary spins CTq and cr|[ play a special role. 
o~% = (Text i s fixed to 1 but <7q = of nt is free (but 
conserved by Ha) and labels the TS. The spectrum of the 
ICTF can thus be studied separately for Cg = +1 or erg = 
— 1. One takes A > (ferromagnetic chain) without loss 
of generality. The sector with erg = +1 thus corresponds 
to unfrustrated boundary conditions for the (pseudo-spin) 
chain. On the other hand, choosing erg = — 1 amounts 
to impose at least one Ising domain wall in the system.EH 
In the thermodynamic limit both sectors have the same 
energy per site but they may have a finite difference in the 
total energy (gap). In the paramagnetic phase (A < 1), 
because of the finite spin-spin correlation length £(A), 
the frustration has an exponentially small effect on the 
ground-state energy and the energy difference between 
the two TS is AE ~ exp(-L/£) with £ ~ Ln(l/A) _1 



IV. MIXING THE T x = ±1 SECTORS WITH 
MONOMERS 

The perturbation Hi discussed in the previous section 
is not the only way to remove the topological degeneracy. 
It is well known that in the presence p£ mobile monomers, 
T x = (— l)^^is-jno longer conservedra This property was 
used in Refsfl'U to mix the topological sectors. We will 
show that in our kagome geometry the arrow represen- 
tation of the dimer model allows to compute exactly the 
spectrum of the system when monomers are allowed to 
be created, to hop and to be destroyed along a line wind- 
ing around the cylinder. As we will see this model is 
closely related to that discussed in the previous section: 
A is replaced by A*, monomers will play the role of the 
visons and the ICTF will have periodic and antiperiodic 
boundary conditions instead of open ones. 



A. Hamiltonian 

We relax the parity constraint if-Tf if = 1 on each 
triangle t so that triangles with one or three incoming 
arrows are allowed. When a triangle has one incoming 
arrow, it is naturally interpreted as the presence of a 
monomer (or hole) at the site of this arrow (see Fig. ||). 
When a triangle has three incoming arrows, we interpret 
it as a monomer, and a dimer which are delocalized over 
the three sites £3 Flipping one arrow (i.e acting with rf ) 
on a dimer state therefore creates two monomers on the 
nearby triangles (one of which may be delocalized over 
three sites). We consider the following Hamiltonian: 



Hn 



^£(i 



—X ^x \ 
T+ T+ TV ) 

*0 *1 t 2 1 



(14) 
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FIG. 5: Top: Kagome lattice in the vicinity of the cut A* 
(dashed line). Bottom: Mixed dimer- monomer configuration 
and the arrow representation. 



En- 

h i=l 



hi 



•tf£(i-«*s) as) 



where h\...Q are the sites around hexagon h and £1,2,3 the 
sites of the triangle £ (see Fig. ||). U is a large energy 
enforcing the constraint on low-energy states. This type 
of model was first considered by Kitaev.B For U > the 
ground-state of the Hamiltonian Eq. |l5| is the same as 
for U = 00 (equivalent to Eq. [l]) since r^r^rf commutes 
with all the a x . However, static pairs of monomers are 
present in excited states with energies greater than 2U 
above the grounchatate. 

As Ioffe et cz,L ,Lrllil we wish to use these monomers to 
couple the two TS. For this purpose the monomers are 
allowed to be created and to propagate along one closed 
loop A* winding around the cylinder (Fig. ^|). The sim- 
plest term which does this is 



-IxU 



E 



(16) 



(this choice of normalization will make the analogy with 
the A perturbation clearer). When a t* term acts on 
a site located between two triangles satisfying the con- 
straint {rf rf Tf 2 — 1), it creates a pair of monomers. 
When it acts on a pair of triangles violating the con- 
straint, the pair of monomers is destroyed. If it acts on 
a site located between triangles with different values of 
T t T ti T t- 2 i a monomer hops from one triangle to the other. 



B. Mapping to the ICTF 

As in Eq. ||, H'p^r H[ splits into a bulk and a one- 
dimensional parts H 



H' + H\ = H' hulk + Ha* 



(17) 



Hi 



bulk 



£^ + [/£(l-r t > t X 2 ) (18) 



Ha* 



f EC 1 

teA* 



rtjlrl) (19) 



One can simply check that Ha* and 7Y(j ulk commute with 
each other so that one has to study a one-dimensional 
model Ha*- This model is identical to a closed ICTF 
(with periodic or antiperiodic boundary conditions) as 
explained below. Each triangle £ crossed by A* cor- 
responds to a site of the spin chain. The associated 
transverse-field term for this Ising spin is: 



"t ~ 'to'ti'ta 
We define the z component of the spins as 



'0 '1 



i(t) 



(20) 



(21) 



where is an (arbitrary) origin on A* and i(t) is the 
site of K in common with triangles £ and £ — 1. It is 
simple to check that the a x and a z defined above (not 
to be confused with a x and a z ) obey the usual Pauli 
matrix algebra and play the role of spin-^ operators. In 
addition, these definitions insure that 



a t a. 



t+i 



(22) 



where £3 is the common site between triangles £ and £ + 1 
(as in Fig. ||). Because A* is a closed curve, a special 
care is needed for the last term: 



'L-1 U 



f (23) 

where T z (defined in Eq. |]) commutes with Ha* (in the 
same way we had before [Ha,T x ] =0). It can be suc- 
cessively set to ±1 to obtain the whole spectrum. With 
these notations TIa* reads 

r T HA* = -//(TX-i^ + E^+i) 



U 



L-l 



- E °t + cst - 



(24) 



t=o 



This is the Hamiltonian of an ICTF with periodic bound- 
ary conditions when T z = 1 and antiperiodic bound- 
ary conditions when T z = — 1. As for Ha, the model 
has, in the thermodynamic limit, a phase transition at 
fi c = 1 between a paramagnetic phase and a ferromag- 
netic phase. The calculation of the gap between the two 
sectors amounts to study a closed ICTF with periodic 
and antiperiodic boundary conditions. The exact result 
for the energy difference AE' is derived in the Appendix 
[b| (see also Ref. |lj) with the help of a Jordan- Wigner 
transformation. In the limit of a large system (L 1) it 
is given by 



AE' 



E, - E_ 



~ 2U 



1 - 



HA* 



Lit 

~ 2U(n- 1) 



P [i L for n<l 



for ji > 1 



(25) 
(26) 
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L=10 




FIG. 6: Effect of a change of boundary conditions on the 
ground-state energy of an ICTF with N = 10 spins as a func- 
tion of the (ferromagnetic) exchange A. Circles: closed chain 
with periodic and anti-periodic boundary conditions. Thin 
line: Ground-state energy difference between open chain with 
ferromagnetic (t " ■ T) an d antiferromagnetic (| ■ ■ ■ T) bound- 
ary conditions. Thick line: N — oo case (same result for the 
open and closed chains). 



The result for L = 10 is plotted in Fig. ^. As for Ji a , the 
critical value [i c = 1 separates a regime with an exponen- 
tially small splitting between the sectors (E q. |25| ) and a 
regime with a finite gap between them (Eq. |26|). In the 
spin language, the ferromagnetic phase is characterized 
by (ct z ) 7^ 0. Going back to dimer/monomers variables, 
we find that of flips all the arrows located on A* between 
the origin and t. Thus, it creates or annihilates a pair of 
monomers sitting at both ends of the string or move a 
monomer from one end to the other. In both cases of cre- 
ates or destroys a monomer on the triangle tS3 (<r z ) =/= 
can thus be interpreted as a condensation of monomers 
along A*. This is equivalent to the condensate of visons 
mentioned in the case of TIa when A > 1. 

To conclude this section we discuss a difference be- 
tween the A and the /i perturbations. In the limit where 
A — > oo no dimer can sit on A anymore and the torus 
is reduced to a rectangle, as if the lattice had been cut 
with scissors. If the same geometrical picture was true 
for the perturbation fx along A*, one would erroneously 
conclude that the system is effectively transformed into 
two cylinders with a 2 x 2 = 4-fold ground-state degener- 
acy. This is incorrect for the following reason. The path 
A* chosen to define T Z (A*) (Eq. g) may be shifted by 
multiplication with a a x operator: 



a x {h)T z {A\) — T Z (A* 2 



(27) 



where h is an hexagon next to AJ and A 2 is identical to 
A\ except that it passes on the other side of h. From 
Eq. H?], we see that T z (A^)|0) does not depend on the 
location of A£ as long as this closed curve is deformed 
by passing only on hexagons with a x (h) = 1. In the 
case of a the perturbation /i, the ground-state satisfies 
a x (h)\0) = |0) for all hexagons h, even those along A* 
(this is not true for the A perturbation which does not 
commute with o~ x ). Therefore one cannot independently 



flip the topological sector of the "upper" cylinder with 
some T* without changing the sector of the "lower part" 
which is controlled by a T|, since Tf |0) = T||0). 



V. A TOY MODEL FOR A TOPOLOGICAL 
QUBIT 

Kitaevi suggested that systems with topologically de- 
generate ground-states could be used to realize qubits 
protected from decoherence. Xivis, suggestion was &<fx\ 
made more precise by Ioffe et aZ.LrEil and Dougotef aZJl3'E2l 
who proposed to use Josephson junctions arrays to imple- 
ment such a system. As mentioned in the introduction, 
the topological nature of the degeneracy makes it difficult 
to "manipulate" (perform unitary rotation) because it is 
almost insensitive to local couplings. On the other hand, 
it may be difficult to apply a perturbation correspond- 
ing to a non-local operator (such as T z or T v ). If how- 
ever the system allows for atuhardware implementation 
of such a non-local operator JEl it represents a dangerous 
channel through which perturbations could bypass the 
topological protection and contribute to decoherence of 
the qubit. If such a non-local coupling to the system ex- 
ists, one must be able to "disconnect" it efficiently when 
it is not active. 

The clever solution proposed in Ref. g,|ll] consists in 
perturbing the system by two external potentials acting 
along two lines, exactly as A and A*. While these per- 
turbations are local (more precisely they are sums of local 
terms), they induce a splitting of the two ground-states 
proportional to A L and therefore induce a slow preces- 
sion of the qubit. In this section we take advantage of 
the exact solution of the model to discuss these effects 
beyond the regime where A L <C 1. 



A. Unitary rotations 

Combining the perturbations along A and A* the 
Hamiltonian is 



"Ho' + %v + Hi 

6 

-En^+^Et 1 -^^) 



h i=l 



(28) 



(29) 



ieA' 



ieA 



From the previous calculations we know that TL{\, /i = 0) 
lifts the degeneracy of the two TS. It acts in this two- 
dimensional su bspac e as AE(X)T X , where AE(X), given 
in Eqs. A13 and A15 , is the energy difference between fer- 
romagnetic and antiferromagnetic boundary conditions 
for an ICTF with exchange A (and a unit transverse 
field). On the other hand, H(\ — 0,/i,) mixes the two 
sectors. Its action is described by UAE'(fi)T z where 
AE' (Eq. B23) is energy difference between periodic and 
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antiperiodic boundary conditions foe. an ICTF (with ex- 
change fiU and transverse field 11)^3 

If the system is operated below the critical values of A 
and /z, the qubit precesses at a frequency which is expo- 
nentially small in the system size. The cylinder topology 
and the low density of monomers protect the degener- 
acy of the spectrum. This is the regime mentioned by 
Ioffe et al. However, the system size cannot be too large 
because the time required for a unitary rotation would 
become exponentially long. On the other hand, if the 
system is not large enough, its topology does not fully 
protect it from decoherence by external perturbations. 

If one of the external parameters (A or ju,) is pushed 
above its critical value, the frequency becomes finite, even 
in case of large system size. We may therefore take ad- 
vantage of the phase transitions in a large system. In such 
a case the qubit is topologically protected as long as A 
and n are smaller than their critical values, even if they 
are not precisely set to zero or if they introduce some 
noise. It is only during the "manipulation" (A(f) > 1 
or /j,(t) > 1) that the state of the qubit evolves (and is 
sensitive to the external noise entering through A or A*). 

To preserve an adiabatic evolution one must avoid 
transitions to other eigenstates. However the gap in the 
spectrum of the ICTF becomes small (of the order of 
~ 1/L) in the vicinity of the transition. This limits the 
typical time of the unitary rotation to be at least of the 
order of L. This linear dependence in the system size is 
an interesting property because it should be compared 
to the exponential dependence (~ A~ L ) present in the 
perturbative regime. Also because of this small gap, the 
temperature must be T <C 1/L to avoid thermal excita- 
tions when A (or fi) is close to 1. Using the transition 
to perform unitary rotation therefore seems to improve 
the time of operation and could enable to use a larger 
system and to benefit from a stronger topological pro- 
tection. It also requires to work at lower temperature 
than for a qubit operated in the perturbative (A,/i -C 1) 
regime only and this may represent a severe limitation. 



B. Reading out the state of the qubit 

We assume that the qubit is in a linear combination 
of the two topological sectors: = a\+) + where 
T x \+) = |+) and T x \—) = — |— ). We wish to measure 
\a\ 2 with a local observable. This is not directly possible 
if the system is very large since any local observable has 
expectation values in |+) and |— ) which are exponentially 
close. Likewise, a local observable has a vanishing matrix 
element between |+) and |— ). A possible procedure could 
be to switch adiabatically the potential A above the tran- 
sition. For a strong enough A, the state |+) evolves to a 
superposition of dimer coverings with no dimer crossing 
A. On the other hand, |— ) evolves to a superposition 
of dimer coverings with one dimer crossing A. This is 
because the parity (T x ) is a conserved quantity under 
the evolution but the ground-state has to minimize 



as X(t) grows. A (projective) measurement detecting the 
presence of a dimer on some bond crossing A will thus 
give 1 with probability \a\ 2 /L. The whole operation has 
to be repeated a large number of times (~ L) before \a\ 2 
is known with a reasonable accuracy but one may im- 
prove the efficiency of the measurement by having a bond 
along A where the energy cost of a dimer is less than on 
other bonds (in which case the dimer, if present, will lo- 
calize on this particular bond when A becomes large). 
Of course the reading procedure described above suffers 
from the same limitations (time proportional to L and 
low temperature) as the unitary rotation. 



VI. CONCLUSIONS 

We have shown that the QDM of Eq. [I] can be simply 
solved in the presence of two kinds of perturbations: an 
external potential that couples to dimers crossing a line 
or the inclusion of monomers. This provides a simple ex- 
ample of system with a Z 2 fractionalized phase where the 
topological degeneracy is destroyed by tuning an external 
parameter through a quantum phase transition (belong- 
ing to the classical Ising 2D universality class). 

We also discussed some properties of this toy model 
from the point of view of an ideal topological qubit, in 
which case the exact solution allows to follow the two low- 
est eigenstates as a function of some external parameters. 
These two parameters can be used to perform unitary ro- 
tations of the qubit and provide an exactly, solvable ver- 
sion of some ideas introduced previously.E In addition, 
we pointed out that the phase transition could, in prin- 
ciple, be used to improve the robustness to decoherence 
because it could enable to use a larger (although not infi- 
nite) system. Concerning the measurement of the qubit, 
we also emphasize the interesting properties of the phase 
transition as it turns a non-local property (parity of the 
number of dimers crossing a line) into a local property 
(dimer density). From this point of view, we note that 
the method has some resemblance with the flux trapping 
experiment imagined by Senthil and FishcrtL2l to detect 
visons in a Z 2 fractionalized systems. 
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APPENDIX A: GROUND-STATE ENERGY OF 
AN ICTF WITH OPEN BOUNDARY 
CONDITIONS 



The third equation (A10) comes from the boundary at 
ri = L + 1 and is a constraint on the available us. 



We apply fixed boundary conditions to an open ICTF 
and compute the energy difference between the case of 
ferromagnetic boundary conditions (two fixed up spins 
at the ends) and antiferromagnetic boundary conditions 
(one up spin at one end and a down spin at the otheci. 
This result has been obtained recently by Dougot et alS3 
but we give here for completeness a detailed derivation 
of the result. 



Hamiltonian and free fermions 



The Hamiltonian is 



n 



(Al) 



n=0 



with two fixed spins at the ends of the chain: er£ +1 = 1 
and (Tq = ±1 depending on the boundary conditions. 
Using a standard Jordan- Wigner transformation to rep- 
resent the Ising operators with spinless fermions: 



u n — 1c n Cn 



n-l 



2<4i exp iir ^ etc. 



i=o 



'n u n+l 



(4 +Cn)(Cn+l ~4+l) 



(A2) 
(A3) 
(A4) 



TL is quadratic in the fermion operators and can be di- 
agonalized by a Bogoliubov transformation. To find the 
quasi-particle creatioit-operators S we consider the fol- 
lowing form ( Ansatz) :t-l 



(A5) 



L + l L + l 

ft = £ W »(4+c„)+t&H5>»(4 -c„)(A6) 



n=0 



n=0 



where u> and b(u>) have to be determined. One can check 
that 



[n.di] = e(u)4 

with E(lo) > 

provided the following equations are satisfied: 

E(u) = -2ib{uj){l + Xuj- 1 ) 
iE(uj)b(uj) = -2(1 + Aw) 



(A7) 



.L + l 



-L-l 



-A(« 



,L+2 



-L-2\ 



(A8) 
(A9) 
(A10) 



The t wo fi rst equations come from the terms < n < L 
in Eq. A7 . These equations determine the energy of the 
quasi-particles as a function of u: 



E(u) = 2VA 2 + A(w + w- 1 ) + l 



(All) 



2. Case A > 



For A > j±i the Eq. |A10| has L + l distinct solutions 



of the form: 



A 



e z " with fce]0,?r] 

sin ((L + l)fc) 
~sin((L + 2)fc) 



(A12) 



These fermionic levels are those required to describe the 
L + l spins n = {)■■■ L. Because we have chosen E to be 
always positive, the absolute ground-state (whatever (Tq) 
is the vacuum |0) of the d'^. The lowest excited state is 
d' ko |0) where we have added the fermion with the smallest 
energy A. It corresponds to the solution fco of Eq. A12 
which is the closest to tt. This solution can be calculated 
by an expansion in 1/L and we obtain 



2(A-1) + 



Att 2 



(A-l)L 2 



+ 0{L- i ) (A13) 



So far we have not specified the value of <Tg correspond- 
ing to each level. In the limit where A 3> 1 it is clear 
that the ground-state is in the sector <7q = 1. Since one 
can show that no level-crossing occur for A > in such a 
finite chain, the fermion vacuum |0) satisfies erg |0) = |0) 
for all A > and correspond to a state of the system with 
ferromagnetic boundaries. On the other hand, inserting 
any d^ fermion changes the sign of ctq (since all the dj, 
anticommute with Ctq = cj, + Co). The first excited state 
d k |0) is thus the ground-state of the system with anti- 
ferromagnetic boundary conditions (fig = —1) and the 
gap between the two sectors is given by A (Eq. A13). 



3. Case < A < £±| 



In the range > A > jrQ, only L solutions of the 



form Eqs. A12 exist. The "missing" solution uiq has the 
lowest energy and is real: lo e] — oo, 0]. It corresponds to 
a bound-state (imaginary wave- vector) for the fermions. 
In the thermodynamic limit one has ujq = — \ and finite- 
size corrections can be evaluated: 



LOq 



+ A 



2L+1 



1-A 2 



0(LX 



4L\ 



which gives an energy gap 



2A 



L+l 



1 — A 



0(L\ 



3L\ 



(A14) 



(A15) 



As before, this gap is the ground-state energy difference 
between antiferromagnetic and ferromagnetic boundary 
conditions for the spin chain. 
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APPENDIX B: GROUND-STATE ENERGY OF 
AN ICTF WITH PERIODIC OR ANTIPERIODIC 
BOUNDARY CONDITIONS 

We compute the ground-state energy of a closed ICTF 
with periodic and anti-periodic boundary conditions. 
The latter result has also been obtained recently by 
Doucot et alt3 but we give here a detailed derivation 
of the result. We also note that this calculation has 
some similarities with the evaluation of the ground-state 
energy splitting in the triangular lattice QDM at r-the 
Rokhsar Kivelson point (using a Pfaffian technique) a 



Using the explicit form of the transformation and the 
anti-periodic boundary conditions on the fermions, one 
can show that the vacuum |0) of the d\ fermions satisfies: 



cxp (iir^t&Cn J|0) = +|0) 

V n=0 / 



(B8) 



This is consistent with Eq. B4 and the ground state is 
thus |0). Its energy is 



E F 



(B9) 



1. Periodic boundary conditions 



The Hamiltonian is : 

L-l 

n = - J2 < 



L-l 

V E a n a n+l 
n=0 



(Bl) 



with cr£ = ctq. The Ising op erators are r epr esented with 



spinless fermions as in Eqs. A2, A3 and A.4. Due to the 
periodic boundary conditions we have also 



^L-i^o = -( c l-i + c L -i)(c - 4) 
x exp in ^ 4i c n 



n=0 



It is simple to check that 

L-l / L-l 

II °n = ex P i7T C * Ci 



(B2) 



(B3) 



71=0 



n=0 



is a conserved quantity. The spectrum can thus be stud- 
ied separately in the sectors rin=o °n = il- However H 
has off-diagonal matrix element in the natural Ising basis 
which are all < 0. This insures that the ground-state has 
only positive weight in this basis and it therefore belongs 
to sector 



L-l 



i 



(B4) 



In the following we thus consider fermions subjected to 
anti-periodic boundary conditions (se e Eqs. [B3| and B4 
and the — sign in the r.h.s of Eq. B2). 

After Fourier transform the Hamiltonian becomes: 



K 



ch am 



[(«Msin(fe)4c^ fe +H.c 



-2c\c k (/ucos(fc) + 1) + 1 



(B5) 



which is diagonalized by a Bogoliubov transformation 

1 



n, 



ch a 



d\d k 



s(k) = 2v//i 2 + 1 + 2ficos(k) 



(B6) 
(B7) 



2. Anti-periodic boundary conditions 

To insure that ct£ — — Oq, the fermions are now sub- 
jected to periodic boundary conditions (see Eq. [B2] ). 
However, for /i > 1 it is necessary to add one dt 
fermion to insure the correct parity under a global spin 
flip. Since the dispersion relation e(k) is minimum in 
e(k = 0) = 2\fi — 1| the ground-state for fi > 1 is 



li) = 4lo) 



(BIO) 



The ground-state energy of the chain with anti-periodic 
boundary conditions is thus 

E A = -- e ( fc ) + 2 (a* - 1) for /i > 1 (BH) 

= -- ^ e(fc) for 0<^<1 (B12) 



3. Energy difference 

From the calculation above the energy difference be- 
tween the ground-states of the two boundary conditions 
is (for fx < 1) : 

L-l 

E A - E P = ^2 [^/coshao - cos(fc„ + i ) 



n=0 



•\/coshao — cos(fc„) 



where 



and da is defined by 



2nir 



cosh ao 
eta 



M 2 + l 
2fi 

-Mm) 



(> i) 



(B13) 



(B14) 



(B15) 
(B16) 
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Re[z] 



FIG. 7: Contour in the complex plane used to define the 
integral Ic- The crosses on the real axis indicate the poles of 
1/ sin(Lz) and the fat segments on the imaginary axis indicate 
the branch cuts of f(z). The region B is sent to lm[z] = ±oo. 



The difference between the two sums can be related to a 
contour integral Iq in the complex plane: 



where 



and 



E A - E P = y% L h 



h = - 



i 



2in J c sin (Lz) 



dz 



f(z) = ^coshevo — cos(z) 



(B17) 



(B18) 



(B19) 



The contour is shown Fig. [?]. The equality Eq. B17 



be demonstrated by using the fact that the poles inside 



the contour are located at z = k n and z 



and 



have alternating residues proportional to — f(k n )/L and 
+f(k n+ i)/L. The contour can be decomposed into sev- 
eral regions: I Q = Ia + Ia' +Ib+Ic + Ic (Fig. 0). Using 
the odd parity of the integrand under z —* — z one has 



Ic = Ic and using the periodicity under z — > z + 2w one 
finds I a + I A' — 0- The integrand is exponentially small 
when lm[z] — > ±oo so the region B does not contribute. 
We therefore have Iq = 21 c- The integral over the region 
C is given by the discontinuity of the integrand along the 
branch cut: 



1 

2iTT 



idr 



f{ir + o ) 
sin(L(ir + o - )) 



f(ir + o+) 



sin(L(ir + o+)) 



^VcosMO-coshM 
sinh(Lr) 



When the system size is large (L — > oo) the behavior of 
Ic is dominated by values of r close to ccq. In this limit 



V sinh an 
Vsmhao 



dre r \Jr — an 



-Lu a r - 



L~2 



(B21) 



so that the energy difference is 

E A -Ep = 2^2^ Lie 



2u sinh an r „ 

2\l— ^ e - Lao 

Lit 



Lit 



ii 2 



(B22) 
(B23) 



The calculation of Ea — Ep for fi > 1 is almost iden- 
tical to the n < 1 case described above. The difference 
between the sums of e(k) on 'even' and 'odd' momenta is 
again expressed with the i nteg ral I q bu t with ao = m(/i). 
Combining this with Eqs. B9 and Bll gives: 



E A -E P = 2(/i - 1) + 2 



1 



Ltt 



— L 



(B24) 
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In fact H can be any trivalent lattice (each site has three 
neighbors) with any topology: disk, sphere, cylinder, torus 
etc... As an example, H can be an hexagonal lattice with 
periodic boundary conditions in one direction (cylinder). 
Then one constructs a new lattice K by the medial lattice 
construction: sites of K are the middle points of the bonds 
of H .E-3 In the following, for simplicity, we will choose H to 
be an hexagonal lattice and K as a kagome lattice (hence 
the names H and K). 

In terms of dimers, is the sum of all possible "kinetic 
energy" terms (ring exchange) which can be defined inside 
the star surrounding h. Those terms allow from 3 to 6 
dimers to move. n 
Those bond variables are Z2 gauge degrees of freedomu 
while the leeal constraint is the Gauss law of the gauge 
description .113 As for <r^ = Ylt=i T i > ^ * s * ne 6 au 6 e nux 
going through the hexagon h. 

They are independent pseudo-spin operators. Notice how- 
ever that if the system has no edge (torus or sphere), they 
are subjected to the an additional constraint W h o~ x (h) ~ 
1. This comes from the fact that, in such a geometry, 
Y\ h <y x {h) flips every arrow twice and thus reduces to the 



identity. 

23 Here local means that the associated loop does not wind 
around the whole system (cylinder). 

24 The case where the potential term pins the dimers along 
some reference configuration in the whole system is dis- 
cussed in Ref. ^ and is equivalent to a 2D Ising model in 
transverse field. 

25 These are bonds where the arrows have different orienta- 
tions in c and in the reference. 

26 The relation between the r and a operators is the same as 
the standard duality between Ising models and Z2 gauge 
theory in 2+1 dimensions. See Ref. ^ 

27 The relation between topological sectors of QDM and 
boun|d|SBj>^-ponditions in spin models was already noted in 
Refs.HOO, together with the connexion between confine- 
ment transition and Ising transition. 

28 If a monomer winds around the cylinder it connects a con- 
figuration T x = 1 with a configuration T x — — 1. 

29 In the 2,2 gauge theory language the constraint is the Gauss 
law and allowing triangles with t^t^t^ = — 1 amounts to 
allow gauge charges ( "matter" ) in the system. 

30 U needs not have the same value in the bulk and along A*. 
One can also choose U = 00 in the bulk and U < 00 along 
A* so that monomers can exist and propagate only along 
the the chain. 

31 This should be compared with a\ which creates a vison on 
hexagon h. 

32 In the circuit proposed in Ref. ul the state of the qubit 
could be measured through a weaK Josephson junction con- 
necting the tow edges of the cylinder. 

33 Since any unitary operation can be performed using T x and 
T z sequentially, there is no need to analyze the spectrum 
with A / and /j, 7^ simultaneously. 

34 Eq. will insure that <tq = cj + Co do not appear in 
so that dlj and anticommute with oq. 



